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ABSTRACT 

Recent Chandra observations of clusters of galaxies revealed the existence of 
a sharp ridge in the X-ray surface brightness where the temperature drops across 
the front. This front is called the cold front. We present the results of two- 
dimensional magnetohydrodynamic simulations of the time evolution of a dense 
subcluster plasma moving in a cluster of galaxies. Anisotropic heat conduction 
along the magnetic field lines is included. In the models without magnetic fields, 
the numerical results indicate that the heat conduction from the hot ambient 
plasma heats the cold dense plasma of the subcluster and diffuses out the cold 
front. When magnetic fields exist in a cluster of galaxies, however, cold fronts 
can be maintained because the heat conduction across the magnetic field lines 
is suppressed. We found that, even when the magnetic fields in a cluster of 
galaxies are disordered, heat conduction across the front is restricted because the 
magnetic field lines are stretched along the front. Numerical results reproduced 
the X-ray intensity distribution observed in the A3667 cluster of galaxies. 

Subject headings: conduction — MHD — galaxies: cluster: magnetic fields — X-ray: 
galaxies — plasmas 



1. INTRODUCTION 



High-resolution Chandra observations have revealed various substructures in the hot 
X-ray emitting plasmas in clusters of galaxies. For example, sharp discontinuities of X- 
ray surface brightness distributions have been found around subclusters moving in a cluster 
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of galaxies such as A2142 (Markevitch et al. 2000) and A3667 (Vikhlinin, Markevitch, & 
Murray 2001b). The temperature drops sharply across this boundary between the grav- 
itationally confined, cold dense plasmas in the moving subcluster and the hot, less-dense 
ambient plasma. This sharp discontinuity is called the cold front. It is not a shock front, 
because Chandra images of the X-ray surface brightness show that the temperature decreases 
in the denser side. 

Ettori & Fabian (2000) suggested that the sharp temperature gradient observed by 
Chandra in A2142 requires the heat conductivity near the front to be reduced by several 
orders of magnitude from the classical Spitzer value (Spitzer 1962), ^sp = 5 x 10~^T^/^ 
ergs~^ cm~^ K~^. 

Based on the observations of A3667, Vikhlinin, Markevitch, & Murray (2001a) suggested 
that turbulent, intergalactic magnetic fields stretched along the cold front reduce the growth 
rate of the Kelvin- Helmholtz (K-H) instability, which possibly destructs the front. Further- 
more, heat conduction across the magnetic fields can be suppressed. It has been recognized 
that heat conduction in a cluster of galaxies can be very efficient (e.g., Takahara & Ikeuchi 
1977) because the Coulomb mean free path, Zc ~ 4 (A;T/5 keV)^(n/10~^cm~^)~^ kpc, is large. 
Here, n is the number density. The time required for heat to diffuse by conduction across 
a length L is given by rgp ~ pL'^ /nsp ~ 10^(A;T/5 kcV)-^'^/2(n/10--'^ cm-3)(L/100 kpc)^ yr, 
where p is the density. However, if magnetic fields exist, the characteristic length of the heat 
exchange across magnetic field lines is reduced to the Larmor radius, Tl = 2500 {B/1 /iG)^^ 
(T/5 keV)^/^ km. Since the typical magnetic field strength in a cluster of galaxies is ~ pG 
(e.g., Kronberg 1994; Carilli & Taylor 2002), the heat conductivity across magnetic field 
lines is drastically reduced. 

In this Letter, we present the results of numerical simulations of the formation of cold 
fronts. These simulations include the effects of heat conduction. We adapted the simulation 
code, originally developed for solar ffare simulations (Yokoyama & Shibata 1997, 2001), 
which includes anisotropic heat conduction along magnetic field lines. 



2. SIMULATION MODELS 

We simulated the time evolution of the cluster plasma in a frame comoving with the 
subcluster. The basic equations are as follows, 

|^ + V.(pt;) = 0, (1) 
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where p, v, p, B, and ip are the density, velocity, pressure, magnetic fields, and gravitational 
potential, respectively. We use the specific heat ratio 7 = 5/3. The subscript || denotes the 
components parallel to the magnetic field lines. We assume that heat is conducted only along 
the field lines. We solved the magnetohydrodynamic (MHD) equations by a modified Lax- 
Wcndroff method with artificial viscosity. The heat conduction term in the energy equation 
is solved by the SOR method (see Yokoyama & Shibata 2001 for details). The units of 
length, velocity, density, pressure, temperature, and time in our simulations are ro = 250 kpc, 
vq = lOOOkms"^ po = 5 X 10~^'''gcm~^, po = 3 x 10~^^ erg cm"^, Tq = 4keV, and to = 
fo/vo — 3 X 10^ yr, respectively. We adopted two-dimensional Cartesian coordinates {x,y) 
by assuming translational symmetry in the z-direction. The typical size of the simulation 
box and the number of grid points are 2.5 Mpc x 2.5 Mpc, and {Nr^,Ny) = (1201,1201), 
respectively. 

Figure 1 shows the initial density distribution. We assume that the subcluster is a sphere 
having a radius r = 500 kpc, and it is isothermal with temperature kTi^ — 4keV. We assume 
that the gravitational potential at r = (x^ + y^)^/^ is given by tlj{r) = (Cg/7) ln[l + (r/rc)^], 
where we use the core radius Vc = 290 kpc, and the sound speed inside the subcluster 
Csin = 1000 km s^^, where the subscript in denotes the values inside the subcluster. The 
density distribution is determined by the hydrostatic equation Vp = —p'Vip, as Pm{r) = 
Pc/[1 + (r/rc)^], where we adopt pc = 10~^^gcm~^. We also assume that the subcluster is 
surrounded by hot (/cTout = 8keV), less-dense (pout = 0.25 po) plasma, where the subscript 
out denotes the values outside the subcluster. The gravitational potential is taken to be 
fixed throughout the simulation. 

In order to simulate the moving subcluster in a cluster of galaxies, we assume that 
ambient plasmas initially have a uniform speed with Mach number M — Vx/cgout — 1, where 
Csout is the sound speed of the ambient plasmas. The Mach number with respect to the 
sound speed inside the subcluster is M' — Vx/cgin — V^- 
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Table 1 shows the model parameters. Models H and HC are without magnetic fields. 
In models M and MC, we assume uniform magnetic fields parallel to the y-axis. In these 
models, the initial ratio of the ambient gas pressure to the magnetic pressure {(5 = Pgas/Pmag) 
is taken to be /3o = 100, which corresponds to the magnetic field strength 5 ~ 3 /iG. 

In models HC, MC, and MTC we included heat conduction. When the magnetic field 
does not exist, heat is conducted isotropically. When the Coulomb mean free path becomes 
comparable to or greater than the temperature scale height It = T/|VT|, the heat flux Q is 
saturated as Qsat = -[^t/(4.2 /c)]/tspVT when It < 4.2 /c (e.g., Cowie & McKce 1977). In 
our initial model, since Ic — lOkpc and It < 6kpc around the cold front, Qsat at the front is 
smaller than 15% of the classical value. However, this reduced heat flux is large enough to 
smear the cold front within 10^ yr because the conduction time scale across the cold front is 
T« ~ {A.21JIt)pIt/ks^ = 2 X 10^ (A;r/5keV)-5/2(n/10-=^cm-3)(/5./6kpc)(4.2/e/42kpc)yr. 

When the magnetic field exists, heat is conducted only in the direction parallel to 
the magnetic field lines. The conductivity parallel to the field fines is k\\ — b x 10~^T^/^ 
ergs~^ cm~^ K~^. Conversely, the conductivity across the field fines is k±_ — 0. The con- 
duction time scale along the field lines is ~ pL"^ / ~ 10''' (/cT/S keV)~^/^(n/10~^ cm~^) 
(L/100kpc)Vr- 

In model MTC, we simulated a subcluster moving in an ambient plasma with initially 
disordered magnetic fields, and we included the effects of heat conduction. In this model, we 
superposed magnetic fields created by thin currents fiowing in the 2;-direction. The location 
and strength of the thin currents are chosen randomly. In this particular model, the size of the 
simulation box is 7.5 Mpc x 2.5 Mpc, and the number of grid points is {N^, Ny) — (901, 301). 

The left boundary at x = — 5 is taken to be a fixed boundary, except in model MTC. In 
model MTC, the magnetic fields at the left boundary are extrapolated from magnetic fields 
in internal points. Other quantities are fixed at the left boundary. The effects of the left 
boundary condition is negligible in the time range we report in this Letter. Other boundaries 
are free boundaries where waves can be transmitted. 

3. RESULTS 

3.1. Effects of Magnetic Fields on Heat Conduction 

Figure 2 shows the distributions of density, temperature and X-ray intensity for model 
MC (upper panels) and model HC (lower panels) at t = 3.0 (0.9 Gyr). Model MC includes 
the effects of anisotropic heat conduction; that is, heat is conducted only in the direction 
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parallel to the magnetic field lines. The X-ray intensity distribution is visualized from the 
simulation results as the logarithm of p^. In model HC with isotropic heat conduction, the 
subcluster plasma is heated by the ambient hot plasma and evaporates rapidly, because heat 
is conducted from the ambient plasma to the cold subcluster plasma. Thus, the sharp front 
of the temperature distribution disappears. We also carried out simulations that included 
the saturation effects of heat flux, and we confirmed that the effects are negligible. In model 
MC, the heat conduction and plasma diffusion across the front are suppressed because the 
magnetic field lines are parallel to the cold front. The X-ray intensity distribution in model 
MC shows a steeper gradient at the leading edge of the subcluster than that in model HC. 
The numerically obtained X-ray intensity distribution of model MC reproduces the overall 
structures of X-ray intensity observed in A3667. At t = 3.0, the magnetic field strength near 
the front increases by a factor of 8, and j3 decreases by a factor of 100 from the initial state. 
In this model, the lifetime of the cold front is more than 10^ yr. 

Numerical results also reproduced the bow shock upstream of the cold front. The 
interaction between the inflowing hot ambient plasma and the cold subcluster plasma drives 
another shock propagating inside the subcluster. This shock is located at {x^y) = (2.5,0) in 
the upper panels of Figure 2. The shock propagating inside the subcluster initially smears 
the cold front by shock heating. Subsequently, adiabatic expansion following the shock 
compression cools the subcluster plasma. This cold plasma is pulled back by gravity and the 
cold front is recovered. 



3.2. Comparison with Models without Heat Conduction 

Figure 3 shows the temperature distribution at t = 3.0 in models M (left) and H (right). 
In these models, heat conduction is not included. By comparing model H with HC, we 
confirm that when the magnetic field does not exist, cold fronts can be maintained only if 
the heat conductivity is small enough. 

The results of model H indicate that the K-H instability grows at the contact discontinu- 
ity between the cold, dense subcluster plasma and the hot, less-dense ambient plasma. The 
growth of the K-H instability is not prominent near the stagnation point {x,y) ~ ( — 1,0), 
but the K-H vortices grow around (x, y) ~ (2, ±2). The K-H instability does not disrupt the 
cold front, because the velocity difference inside and outside of the contact surface is small 
in the upstream side of the dense subcluster. In model M, the magnetic tension force along 
the front suppresses the K-H instability, even in the sides of the dense subcluster, where the 
velocity shear is large. Except for this stability for the K-H instability, the numerical results 
of model M are similar to these of model H. 
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3.3. Effects of Disordered Magnetic Fields 

Figure 4 shows the results of model MTC, in which we initially assumed disordered 
magnetic fields. The top and bottom panels show the time evolution of the density and 
temperature distributions at t = and t = 3.0 (0.9 Gyr), respectively. The white curves and 
arrows show the magnetic field lines and velocity vectors, respectively. Numerical results 
show that the tangled magnetic field lines in the vicinity of the front are compressed and 
stretched along the front. The temperature distributions show that heat is conducted along 
the tangled field lines. As a result, the low-temperature subcluster plasma is heated and 
diffuses along the field lines. Since the tangled magnetic field lines are stretched along the 
front, suppression of heat conduction across the field lines enables the existence of cold fronts. 

4. DISCUSSION 

Heinz et al. (2003) carried out hydrodynamical simulations of the interaction of cold 
subcluster plasma and hot ambient matter. They suggested that a cold front can be formed 
after the passage of shock waves through the subcluster plasma, when the cold plasma, 
initially pushed forward by the shock, falls back. However, in their simulations heat conduc- 
tion and magnetic fields were neglected. In this Letter, we presented the results of numerical 
simulations which included both magnetic fields and heat conduction. We confirmed that, 
since the Spitzer conductivity is very large in hot ambient plasma, cold fronts cannot be 
maintained longer than 0.3 Gyr, unless the heat conduction across the front is suppressed by 
the magnetic fields. 

We also carried out numerical simulations which included heat conduction starting from 
the same initial condition as Heinz et al. (2003). We found that, when we adopt the Spitzer 
conductivity, the subcluster evaporates rapidly and the cold front is not formed. When we 
included magnetic fields and the effects of anisotropic heat conduction, the numerical results 
resembled Heinz's results. As sTiggcsted by Heinz et al. (2003), after the passage of shock 
waves, the cold plasma pulled back by gravity helps sharpen the cold front. 

In §3.2, we presented numerical results that showed the K-H instability does not grow 
near the cold front, because the velocity of inflowing plasmas is reduced after the plasmas 

pass through the bow shock. In the magnetohydrodynamical models, the magnetic fields, 
which stretched along the contact surface, also stabilized the K-H instability. Thus, the K-H 
instability does not threaten the existence of cold fronts. 

In §3.3 we showed results of a model with initially disordered magnetic fields. Numerical 
results show that the magnetic fields near the front are gradually stretched. Thus, the 
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initially tangled magnetic field lines become parallel to the front and suppress the heat 
conduction across the front. We confirmed the suggestion by Vikhlinin et al. (2001a) that 
the stretching of turbulent intergalactic magnetic fields enables the maintenance of cold 
fronts. If the scale of magnetic turbulence is smaller than the mean free path, heat can be 
conducted almost isotropically. The heat conductivity in such a turbulent MHD medium is 
as large as 1/5 of the Spitzer value (Narayan & Medvedev 2001) and still large enough to 
diffuse out the cold front. We expect that coherent fluid motions and the stretching of the 
magnetic fields along the cold front increase the coherent length of field lines and make that 
region less turbulent. 

Let us now discuss the three-dimensional (3D) effects. Since magnetic fields compressed 
in front of the cold front can expand in the third direction, the strength of the magnetic fields 
in this region will be reduced in 3D. However, even such weakened magnetic fields suppress 
the heat conduction across them. Thus, we expect that cold fronts can be sustained even in 
3D. We hope to report the results of such 3D MHD simulations in subsequent papers. 

In this Letter, we assumed that the ambient plasma has magnetic fields. Recently, 
Okabe & Hattori (2003) pointed out the possibility that magnetic fields are spontaneously 
generated by the Ramani-Laval instability around the cold front. Such magnetic fields can 
also help to suppress the heat conduction across the front. 

We thank T. Hanawa, S. Miyaji, N. Okabe, M. Hattori, and U. Pujita for their discus- 
sions. The simulation code CANS that we used was developed as a product of the ACT-JST 
project (PL R. Matsumoto). Wc thank T. Yokoyama for code development. Numerical 
computations were carried out by VPP5000 at the Astronomical Data Analysis Center of 
the National Astronomical Observatory Japan (PI. N. Asai, project-ID yna23). 
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Table 1. Simulation models 
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Fig. 1. — The initial condition of models M and MC. The color contour in the left panel 
shows the density. Arrows show the velocity vectors and white lines show the magnetic 
field lines. The right panel shows the distribution of density (solid curve), pressure (dotted 
curve), and temperature (dashed line) along y = 0. The plasma beta in the ambient plasma 
is /3o = 100. 

Fig. 2. — The effects of heat conduction and magnetic fields. The panels show the dis- 
tributions of density, temperature, and X-ray intensity at t = 3.0 from left to right. The 
upper panels show the results for model MC (with magnetic fields and heat conduction), 
the lower panels show the results for model HC (without magnetic fields and with heat con- 
duction). The X-ray intensity distribution is visualized from the simulation results as the 
logarithm of p^. The white curves and arrows show the magnetic field lines and velocity 
vectors, respectively. 

Fig. 3. — Results of the simulations without heat conduction for models M (left, with mag- 
netic fields) and H (right, without magnetic fields). The color contours show the temperature 
a,t t = 3.0. The solid curves and arrows show the magnetic field lines and velocity vectors, 
respectively. 

Fig. 4. — Time evolutions of model MTC (with disordered magnetic fields) at the central 
region of 5.0 Mpc x 2.5 Mpc at t = and t = 3.0. The top panels show the density distribution 
and the bottom panels show the temperature distribution. The white curves and arrows show 
the magnetic field lines and velocity vectors, respectively. 
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